Sequence Dependence in Nucleosome Dynamics.

The basic packaging unit of eukaryotic chromatin is the nucleosome that contains 145-147 base pair duplex DNA wrapped around an octameric histone protein. While the DNA sequence plays a crucial role in controlling the positioning of the nucleosome, the molecular details behind the interplay between DNA sequence and nucleosome dynamics remain relatively unexplored. This study analyzes this interplay in detail by performing all-atom molecular dynamics simulations of nucleosomes, comparing the human α-satellite palindromic (ASP) and the strong positioning "Widom-601" DNA sequence at time scales of 12 µs. The simulations are performed at salt concentrations 10-20 times higher than physiological salt concentrations to screen the electrostatic interactions and promote unwrapping. These microsecond-long simulations give insight into the molecular-level sequence-dependent events that dictate the pathway of DNA unwrapping. We find that the "ASP" sequence forms a loop around SHL ± 5 for three sets of simulations. Coincident with loop formation is a cooperative increase in contacts with the neighboring N-terminal H2B tail and C-terminal H2A tail and the release of neighboring counterions. We find that the Widom-601 sequence exhibits a strong breathing motion of the nucleic acid ends. Coincident with the breathing motion is the collapse of the full N-terminal H3 tail and formation of an α-helix that interacts with the H3 histone core. We postulate that the dynamics of these histone tails and their modification with post-translational modifications (PTMs) may play a key role in governing this dynamics.

Production of a Monoclonal Antibody for Histone H2b Isoform H2b3b.

H2b3b is one of the histone H2b isoforms that differs from canonical H2b by five to six amino acids. Previously, we identified H3t as the testis-specific histone H3 variant located in histone cluster 3, which is also the site of H2b3b. In this study, we produced monoclonal antibodies against H2b3b, using the iliac rat lymph node method for rat antibody and the immunochamber method for rabbit antibody. Immunoblot analysis confirmed that our antibodies could specifically discriminate between H2b3b and canonical H2b. Moreover, immunostaining revealed colocalization with a testicular stem cell marker, Plzf, but not with a meiotic marker, Sycp. This indicated that H2b3b is expressed in spermatogenic cells before meiosis. Our monoclonal antibodies enable further studies to reveal specific functions of H2b3b during spermatogenesis. We also hope that the established method will lead to the production of antibodies that can identify other H2b isoforms.

Parental histone transfer caught at the replication fork.

In eukaryotes, DNA compacts into chromatin through nucleosomes1,2. Replication of the eukaryotic genome must be coupled to the transmission of the epigenome encoded in the chromatin3,4. Here we report cryo-electron microscopy structures of yeast (Saccharomyces cerevisiae) replisomes associated with the FACT (facilitates chromatin transactions) complex (comprising Spt16 and Pob3) and an evicted histone hexamer. In these structures, FACT is positioned at the front end of the replisome by engaging with the parental DNA duplex to capture the histones through the middle domain and the acidic carboxyl-terminal domain of Spt16. The H2A-H2B dimer chaperoned by the carboxyl-terminal domain of Spt16 is stably tethered to the H3-H4 tetramer, while the vacant H2A-H2B site is occupied by the histone-binding domain of Mcm2. The Mcm2 histone-binding domain wraps around the DNA-binding surface of one H3-H4 dimer and extends across the tetramerization interface of the H3-H4 tetramer to the binding site of Spt16 middle domain before becoming disordered. This arrangement leaves the remaining DNA-binding surface of the other H3-H4 dimer exposed to additional interactions for further processing. The Mcm2 histone-binding domain and its downstream linker region are nested on top of Tof1, relocating the parental histones to the replisome front for transfer to the newly synthesized lagging-strand DNA. Our findings offer crucial structural insights into the mechanism of replication-coupled histone recycling for maintaining epigenetic inheritance.

Cryo-EM structures of RAD51 assembled on nucleosomes containing a DSB site.

RAD51 is the central eukaryotic recombinase required for meiotic recombination and mitotic repair of double-strand DNA breaks (DSBs)1,2. However, the mechanism by which RAD51 functions at DSB sites in chromatin has remained elusive. Here we report the cryo-electron microscopy structures of human RAD51-nucleosome complexes, in which RAD51 forms ring and filament conformations. In the ring forms, the N-terminal lobe domains (NLDs) of RAD51 protomers are aligned on the outside of the RAD51 ring, and directly bind to the nucleosomal DNA. The nucleosomal linker DNA that contains the DSB site is recognized by the L1 and L2 loops-active centres that face the central hole of the RAD51 ring. In the filament form, the nucleosomal DNA is peeled by the RAD51 filament extension, and the NLDs of RAD51 protomers proximal to the nucleosome bind to the remaining nucleosomal DNA and histones. Mutations that affect nucleosome-binding residues of the RAD51 NLD decrease nucleosome binding, but barely affect DNA binding in vitro. Consistently, yeast Rad51 mutants with the corresponding mutations are substantially defective in DNA repair in vivo. These results reveal an unexpected function of the RAD51 NLD, and explain the mechanism by which RAD51 associates with nucleosomes, recognizes DSBs and forms the active filament in chromatin.

Histone variant H2A.Z and linker histone H1 influence chromosome condensation in Saccharomyces cerevisiae.

Chromosome condensation is essential for the fidelity of chromosome segregation during mitosis and meiosis. Condensation is associated both with local changes in nucleosome structure and larger-scale alterations in chromosome topology mediated by the condensin complex. We examined the influence of linker histone H1 and variant histone H2A.Z on chromosome condensation in budding yeast cells. Linker histone H1 has been implicated in local and global compaction of chromatin in multiple eukaryotes, but we observe normal condensation of the rDNA locus in yeast strains lacking H1. However, deletion of the yeast HTZ1 gene, coding for variant histone H2A.Z, causes a significant defect in rDNA condensation. Loss of H2A.Z does not change condensin association with the rDNA locus or significantly affect condensin mRNA levels. Prior studies reported that several phenotypes caused by loss of H2A.Z are suppressed by eliminating Swr1, a key component of the SWR complex that deposits H2A.Z in chromatin. We observe that an htz1Δ swr1Δ strain has near-normal rDNA condensation. Unexpectedly, we find that elimination of the linker histone H1 can also suppress the rDNA condensation defect of htz1Δ strains. Our experiments demonstrate that histone H2A.Z promotes chromosome condensation, in part by counteracting activities of histone H1 and the SWR complex.

Decoding Chromatin Ubiquitylation: A Chemical Biology Perspective.

Since Strahl and Allis proposed the "language of covalent histone modifications", a host of experimental studies have shed light on the different facets of chromatin regulation by epigenetic mechanisms. Initially proposed as a concept for controlling gene transcription, the regulation of deposition and removal of histone post-translational modifications (PTMs), such as acetylation, methylation, and phosphorylation, have been implicated in many chromatin regulation pathways. However, large PTMs such as ubiquitylation challenge research on many levels due to their chemical complexity. In recent years, chemical tools have been developed to generate chromatin in defined ubiquitylation states in vitro. Chemical biology approaches are now used to link specific histone ubiquitylation marks with downstream chromatin regulation events on the molecular level. Here, we want to highlight how chemical biology approaches have empowered the mechanistic study of chromatin ubiquitylation in the context of gene regulation and DNA repair with attention to future challenges.

Regulation of adipogenesis by histone methyltransferases.

Epigenetic regulation is a critical component of lineage determination. Adipogenesis is the process through which uncommitted stem cells or adipogenic precursor cells differentiate into adipocytes, the most abundant cell type of the adipose tissue. Studies examining chromatin modification during adipogenesis have provided further understanding of the molecular blueprint that controls the onset of adipogenic differentiation. Unlike histone acetylation, histone methylation has context dependent effects on the activity of a transcribed region of DNA, with individual or combined marks on different histone residues providing distinct signals for gene expression. Over half of the 42 histone methyltransferases identified in mammalian cells have been investigated in their role during adipogenesis, but across the large body of literature available, there is a lack of clarity over potential correlations or emerging patterns among the different players. In this review, we will summarize important findings from studies published in the past 15 years that have investigated the role of histone methyltransferases during adipogenesis, including both protein arginine methyltransferases (PRMTs) and lysine methyltransferases (KMTs). We further reveal that PRMT1/4/5, H3K4 KMTs (MLL1, MLL3, MLL4, SMYD2 and SET7/9) and H3K27 KMTs (EZH2) all play positive roles during adipogenesis, while PRMT6/7 and H3K9 KMTs (G9a, SUV39H1, SUV39H2, and SETDB1) play negative roles during adipogenesis.

Trimethyllysine Reader Proteins Exhibit Widespread Charge-Agnostic Binding via Different Mechanisms to Cationic and Neutral Ligands.

In the last 40 years, cation-π interactions have become part of the lexicon of noncovalent forces that drive protein binding. Indeed, tetraalkylammoniums are universally bound by aromatic cages in proteins, suggesting that cation-π interactions are a privileged mechanism for binding these ligands. A prominent example is the recognition of histone trimethyllysine (Kme3) by the conserved aromatic cage of reader proteins, dictating gene expression. However, two proteins have recently been suggested as possible exceptions to the conventional understanding of tetraalkylammonium recognition. To broadly interrogate the role of cation-π interactions in protein binding interactions, we report the first large-scale comparative evaluation of reader proteins for a neutral Kme3 isostere, experimental and computational mechanistic studies, and structural analysis. We find unexpected widespread binding of readers to a neutral isostere with the first examples of readers that bind the neutral isostere more tightly than Kme3. We find that no single factor dictates the charge selectivity, demonstrating the challenge of predicting such interactions. Further, readers that bind both cationic and neutral ligands differ in mechanism: binding Kme3 via cation-π interactions and the neutral isostere through the hydrophobic effect in the same aromatic cage. This discovery explains apparently contradictory results in previous studies, challenges traditional understanding of molecular recognition of tetraalkylammoniums by aromatic cages in myriad protein-ligand interactions, and establishes a new framework for selective inhibitor design by exploiting differences in charge dependence.

The role of cryptic ancestral symmetry in histone folding mechanisms across Eukarya and Archaea.

Histones compact and store DNA in both Eukarya and Archaea, forming heterodimers in Eukarya and homodimers in Archaea. Despite this, the folding mechanism of histones across species remains unclear. Our study addresses this gap by investigating 11 types of histone and histone-like proteins across humans, Drosophila, and Archaea through multiscale molecular dynamics (MD) simulations, complemented by NMR and circular dichroism experiments. We confirm and elaborate on the widely applied "folding upon binding" mechanism of histone dimeric proteins and report a new alternative conformation, namely, the inverted non-native dimer, which may be a thermodynamically metastable configuration. Protein sequence analysis indicated that the inverted conformation arises from the hidden ancestral head-tail sequence symmetry underlying all histone proteins, which is congruent with the previously proposed histone evolution hypotheses. Finally, to explore the potential formations of homodimers in Eukarya, we utilized MD-based AWSEM and AI-based AlphaFold-Multimer models to predict their structures and conducted extensive all-atom MD simulations to examine their respective structural stabilities. Our results suggest that eukaryotic histones may also form stable homodimers, whereas their disordered tails bring significant structural asymmetry and tip the balance towards the formation of commonly observed heterotypic dimers.

The 80th Threonine Residue of Histone H3 Is Important for Maintaining HM Silencing in Saccharomyces cerevisiae.

Gene expression in eukaryotic cells is intricately regulated by chromatin structure and various factors, including histone proteins. In Saccharomyces cerevisiae, transcriptionally silenced regions, such as telomeres and homothallic mating (HM) loci, are essential for genome stability and proper cellular function. We firstly observed the defective HM silencing in alanine substitution mutant of 80th threonine residue of histone H3 (H3T80A). To identify which properties in the H3T80 residue are important for the HM silencing, we created several substitution mutants of H3T80 residue by considering the changed states of charge, polarity, and structural similarity. This study reveals that the structural similarity of the 80th position of H3 to the threonine residue, not the polarity and charges, is the most important thing for the transcriptional silencing in the HM loci.

Using NMR diffusion data to validate MD models of disordered proteins: Test case of N-terminal tail of histone H4.

MD simulations can provide uniquely detailed models of intrinsically disordered proteins (IDPs). However, these models need careful experimental validation. The coefficient of translational diffusion Dtr, measurable by pulsed field gradient NMR, offers a potentially useful piece of experimental information related to the compactness of the IDP's conformational ensemble. Here, we investigate, both experimentally and via the MD modeling, the translational diffusion of a 25-residue N-terminal fragment from histone H4 (N-H4). We found that the predicted values of Dtr, as obtained from mean-square displacement of the peptide in the MD simulations, are largely determined by the viscosity of the MD water (which has been reinvestigated as a part of our study). Beyond that, our analysis of the diffusion data indicates that MD simulations of N-H4 in the TIP4P-Ew water give rise to an overly compact conformational ensemble for this peptide. In contrast, TIP4P-D and OPC simulations produce the ensembles that are consistent with the experimental Dtr result. These observations are supported by the analyses of the 15N spin relaxation rates. We also tested a number of empirical methods to predict Dtr based on IDP's coordinates extracted from the MD snapshots. In particular, we show that the popular approach involving the program HYDROPRO can produce misleading results. This happens because HYDROPRO is not intended to predict the diffusion properties of highly flexible biopolymers such as IDPs. Likewise, recent empirical schemes that exploit the relationship between the small-angle x-ray scattering-informed conformational ensembles of IDPs and the respective experimental Dtr values also prove to be problematic. In this sense, the first-principle calculations of Dtr from the MD simulations, such as demonstrated in this work, should provide a useful benchmark for future efforts in this area.

Uncovering New Conformational States of the Substrate Binding Pocket of LSD1 Potential for Inhibitor Design via Funnel Metadynamics.

Lysine-specific demethylase 1 (LSD1) is a promising therapeutic target for cancer therapy. So far, over 80 crystal structures of LSD1 in different complex states have been deposited in the Protein Data Bank, which are valuable resources for performing structure-based drug design. However, among all of the crystal structures of LSD1, the substrate binding pocket, which is the most efficient druggable site for designing LSD1 inhibitors at present, is very similar no matter whether LSD1 is in the apo or any holo forms, which is inconsistent with its versatile demethylase functions. To investigate whether the substrate binding pocket is rigid or exhibits other representative conformations different from the crystal conformations that are feasible for designing new LSD1 inhibitors, we performed funnel metadynamics simulations to study the conformation dynamics of LSD1 in the binding process of two effective LSD1 inhibitors (CC-90011 and 6X0, CC-90011 undergoing clinical trials). Our results showed that the entrance of the substrate binding pocket is very flexible. Two representative entrance conformations of LSD1 counting against binding with the substrate of histone H3 were detected, which may be used for structure-based LSD1 inhibitor design. Besides, alternative optimal binding modes and prebinding modes for both inhibitors were also detected, which depicted that the key interactions changed along with the binding process. Our results should provide great help for LSD1 inhibitor design.

Histone Readers and Their Roles in Cancer.

Histone proteins in eukaryotic cells are subjected to a wide variety of post-translational modifications, which are known to play an important role in the partitioning of the genome into distinctive compartments and domains. One of the major functions of histone modifications is to recruit reader proteins, which recognize the epigenetic marks and transduce the molecular signals in chromatin to downstream effects. Histone readers are defined protein domains with well-organized three-dimensional structures. In this Chapter, we will outline major histone readers, delineate their biochemical and structural features in histone recognition, and describe how dysregulation of histone readout leads to human cancer.

Polycomb Repressive Complex 2 in Oncology.

Dynamic regulation of the chromatin state by Polycomb Repressive Complex 2 (PRC2) provides an important mean for epigenetic gene control that can profoundly influence normal development and cell lineage specification. PRC2 and PRC2-induced methylation of histone H3 lysine 27 (H3K27) are critically involved in a wide range of DNA-templated processes, which at least include transcriptional repression and gene imprinting, organization of three-dimensional chromatin structure, DNA replication and DNA damage response and repair. PRC2-based genome regulation often goes wrong in diseases, notably cancer. This chapter discusses about different modes-of-action through which PRC2 and EZH2, a catalytic subunit of PRC2, mediate (epi)genomic and transcriptomic regulation. We will also discuss about how alteration or mutation of the PRC2 core or axillary component promotes oncogenesis, how post-translational modification regulates functionality of EZH2 and PRC2, and how PRC2 and other epigenetic pathways crosstalk. Lastly, we will briefly touch on advances in targeting EZH2 and PRC2 dependence as cancer therapeutics.

The CRTC-1 transcriptional domain is required for COMPASS complex-mediated longevity in C. elegans.

Loss of function during aging is accompanied by transcriptional drift, altering gene expression and contributing to a variety of age-related diseases. CREB-regulated transcriptional coactivators (CRTCs) have emerged as key regulators of gene expression that might be targeted to promote longevity. Here we define the role of the Caenorhabditis elegans CRTC-1 in the epigenetic regulation of longevity. Endogenous CRTC-1 binds chromatin factors, including components of the COMPASS complex, which trimethylates lysine 4 on histone H3 (H3K4me3). CRISPR editing of endogenous CRTC-1 reveals that the CREB-binding domain in neurons is specifically required for H3K4me3-dependent longevity. However, this effect is independent of CREB but instead acts via the transcription factor AP-1. Strikingly, CRTC-1 also mediates global histone acetylation levels, and this acetylation is essential for H3K4me3-dependent longevity. Indeed, overexpression of an acetyltransferase enzyme is sufficient to promote longevity in wild-type worms. CRTCs, therefore, link energetics to longevity by critically fine-tuning histone acetylation and methylation to promote healthy aging.

Conformational and Interaction Landscape of Histone H4 Tails in Nucleosomes Probed by Paramagnetic NMR Spectroscopy.

The fundamental repeat unit of chromatin, the nucleosome, consists of approximately 147 base pairs of double-stranded DNA and a histone protein octamer containing two copies each of histones H2A, H2B, H3, and H4. Each histone possesses a dynamically disordered N-terminal tail domain, and it is well-established that the tails of histones H3 and H4 play key roles in chromatin compaction and regulation. Here we investigate the conformational ensemble and interactions of the H4 tail in nucleosomes by means of solution NMR measurements of paramagnetic relaxation enhancements (PREs) in recombinant samples reconstituted with 15N-enriched H4 and nitroxide spin-label tagged H3. The experimental PREs, which report on the proximities of individual H4 tail residues to the different H3 spin-label sites, are interpreted by using microsecond time-scale molecular dynamics simulations of the nucleosome core particle. Collectively, these data enable improved localization of histone H4 tails in nucleosomes and support the notion that H4 tails engage in a fuzzy complex interaction with nucleosomal DNA.

Stochastic nucleosome disassembly mediated by remodelers and histone fragmentation.

We construct and analyze monomeric and multimeric models of the stochastic disassembly of a single nucleosome. Our monomeric model predicts the time needed for a number of histone-DNA contacts to spontaneously break, leading to dissociation of a non-fragmented histone from DNA. The dissociation process can be facilitated by DNA binding proteins or processing molecular motors that compete with histones for histone-DNA contact sites. Eigenvalue analysis of the corresponding master equation allows us to evaluate histone detachment times under both spontaneous detachment and protein-facilitated processes. We find that competitive DNA binding of remodeling proteins can significantly reduce the typical detachment time but only if these remodelers have DNA-binding affinities comparable to those of histone-DNA contact sites. In the presence of processive motors, the histone detachment rate is shown to be proportional to the product of the histone single-bond dissociation constant and the speed of motor protein procession. Our simple intact-histone model is then extended to allow for multimeric nucleosome kinetics that reveal additional pathways of disassembly. In addition to a dependence of complete disassembly times on subunit-DNA contact energies, we show how histone subunit concentrations in bulk solutions can mediate the disassembly process by rescuing partially disassembled nucleosomes. Moreover, our kinetic model predicts that remodeler binding can also bias certain pathways of nucleosome disassembly, with higher remodeler binding rates favoring intact-histone detachment.

ARES-Kcr: A New Network Model Utilizing Attention Mechanism and Residual Structure for the Prediction of Lysine Crotonylation Sites.

Lysine crotonylation (Kcr), as a significant post-translational modification of protein, exists in the core histones and some non histones of many organisms, and plays a crucial regulatory role in many biological processes such as gene expression, cell development, and disease treatment. Due to the high cost, time-consuming and labor-intensive nature of traditional biological experimental methods, it is necessary to develop efficient, low-cost and accurate calculation methods for identifying crotonylation sites. Therefore, we propose a new network model called ARES-Kcr, which extracts three types of features from different perspectives and integrates convolutional modules, attention mechanisms, and residual modules for feature fusion to improve prediction ability in this paper. Our model performs significantly better than other models on the benchmark dataset, with an average AUC of 92% in the independent test set, demonstrating its excellent predictive ability.

Asymmetric distribution of parental H3K9me3 in S phase silences L1 elements.

In eukaryotes, repetitive DNA sequences are transcriptionally silenced through histone H3 lysine 9 trimethylation (H3K9me3). Loss of silencing of the repeat elements leads to genome instability and human diseases, including cancer and ageing1-3. Although the role of H3K9me3 in the establishment and maintenance of heterochromatin silencing has been extensively studied4-6, the pattern and mechanism that underlie the partitioning of parental H3K9me3 at replicating DNA strands are unknown. Here we report that H3K9me3 is preferentially transferred onto the leading strands of replication forks, which occurs predominantly at long interspersed nuclear element (LINE) retrotransposons (also known as LINE-1s or L1s) that are theoretically transcribed in the head-on direction with replication fork movement. Mechanistically, the human silencing hub (HUSH) complex interacts with the leading-strand DNA polymerase Pol ε and contributes to the asymmetric segregation of H3K9me3. Cells deficient in Pol ε subunits (POLE3 and POLE4) or the HUSH complex (MPP8 and TASOR) show compromised H3K9me3 asymmetry and increased LINE expression. Similar results were obtained in cells expressing a MPP8 mutant defective in H3K9me3 binding and in TASOR mutants with reduced interactions with Pol ε. These results reveal an unexpected mechanism whereby the HUSH complex functions with Pol ε to promote asymmetric H3K9me3 distribution at head-on LINEs to suppress their expression in S phase.

Small-molecule tools for YEATS domain proteins.

Chromatin reader domains are protein folds that bind to post-translational modifications of histones and other chromatin-associated proteins. Compared to other families of reader domains, the discovery that YEATS domains bind to acylated lysines is relatively recent. Four human proteins harbor a YEATS domain, and each is present in protein complexes that regulate chromatin and transcription (ENL, AF9, YEATS2, and YEATS4). Without chemical tools to enable temporally resolved perturbations, it is often unclear how reader domains contribute to protein function. Here, we will discuss recent progress in developing small-molecule tools for YEATS domains and highlight their usefulness for making biological discoveries.